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Abstract 

In this paper we use the symmetry reduction method to obtain invariant solutions of the 
ideal magnetohydrodynamic equations in (3+1) dimensions. These equations are invariant 
under a Galilean-similitude Lie algebra for which the classification by conjugacy classes of 
r-dimensional subalgebras (1 < r < 4) was already known. So we restrict our study to the 
three-dimensional Galilean-similitude subalgebras that give systems composed of ordinary 
differential equations. We present here several examples of these solutions. Some of these 
exact solutions show interesting physical interpretations. 
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1 Introduction 

The ideal magnetohydrodynamic (MHD) equations are the most basic single-fluid model 
with the Maxwell's equations for describing the interactions between magnetic and pres- 
sure forces within an electrically conducting fluid for long spatial scale and low frequency 
phenomena in plasma physics. There are not many exact solutions in the theory of ideal 
magnetized fluid governed by the system of MHD equations. We present here some exact 
analytic solutions of this system. By exact solution we mean one obtained by a reduction 
of the full set of partial differential equations (PDEs) of the original system to one or more 
ordinary differential equations (ODEs) involving fewer independent variables. This can 
be achieved by applying the symmetry reduction method (SRM) based on group theory. 
This method for reducing the number of independent variables in an equation is to require 
that a solution be invariant under some subgroups of the Lie symmetry group of the given 
system; the so-called G- invariant solution (GIS). 

The paper is organized as follows. Section 2 briefly gives a general description of 
the ideal MHD equations. In Section 3 contains a short summary of the SRM and its 
application to the MHD system. In Sections 4 and 5, some examples of GIS are presented 
and also interpreted by a physical point of view for two configurations of the magnetic 
field: B = (Bi,B2,0) and B = (Bi,B2,B3), respectively. Section 6 summarizes the 
obtained results which suggest further applications. 
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The fluid under consideration is assumed to be ideal and perfectly conductive, that is, 
neglecting all the dissipative and dispersive effects, such as viscosity, magnetic resistivity, 
thermal conductivity and Hall effect. The fluid is assumed to be unbounded and the 
medium is non magnetic, so its magnetic permeability /i is taken to be one. Furthermore 
we suppose that the fluid is described by a perfect gas equation of state. Under the above 
assumptions the ideal MHD model is governed by the system of nine PDEs: 

^ + p-1 Vp + p-^B X J) 

dp ,^ , 
— + 7(V- v)p 
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--Vx(vxB) 
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where d/dt = d/dt + {\ ■ V) is the convective derivative. Here p is the mass density, 
p is the pressure, v = (vi,V2,V3) and B = (31,62,63) are the flow velocity and the 
magnetic field, respectively and 7 is the adiabatic exponent. Note that throughout this 
paper we denoted by ei, 62 and the unit vectors in the direction of the x-, y- and 
z-axes, respectively. The current density is given by the Ampere's law 

J = VxB, (2.6) 

for which the displacement current is negligible since the flow is non-relativistic. All the 
functions u = (p, v, B) depend on time t and coordinates x = {x,y,z). The system 
of MHD equations fl2.1l) - fl2.5p is quasilinear, hyperbolic, and is written in the normal 
(Cauchy-Kowalewski) form. In equation (12.21) . there is no extraneous force other than the 
Lorentz force (i.e., the magnetic force) 

F^ = JxB, (2.7) 

which can be broken down into a magnetic pressure force and a magnetic tension force 

F™ = -^V(|B|2) + (B-V)B. (2.8) 

It is noteworthy to recall that the integral form of Faraday's law (12. 4p 

d 
di 



B-ds =0, (2.9) 



s 

where S is an arbitrary surface moving with the fluid, corresponds to the magnetic flux 
conservation law, i.e., the "Alfven's frozen- in" theorem p], and has some fundamental 
consequences. The magnetic field lines are "glued" to the medium, and transported 
entirely by convection and there is no diffusion of B through the conducting media. So, 
in "flux-freezing" regime, the convective derivative of B/p can be written as |2j : 

^ vV = 0. (2.10) 



dt \ p / \ p 
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This says that B/p evolves in the same manner as the separation 61 between two points 
in the fluid. Notice that the vorticity a; of the flow and B are physically related by: 

(2.11) 

In hydrodynamics where B = 0, if at one instance lo = 0, the flow will be irrotational for 
all subsequent times. So, the vorticity field obeys a transport equation similar to 02.41) : 
the vortex lines are frozen into the fluid matter. However this result does not hold for 
ideal MHD because the term Vx [Fm/p] is in general not equal to zero, and may induce 
vorticity even though the flow is initially irrotational. 
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3 The symmetry of the magnetohydrodynamic equations 



The language of group theoretical methods for studying systems of PDEs is a very useful 
and suitable tool for investigating the main features of numerous problems in various 
branches of mathematical physics. The task of finding an increasing number of solutions of 
systems of PDEs is related to the group properties of these differential equations. Its main 
advantages appear when group analysis makes it possible to construct regular algorithms 
for finding certain classes of solutions without referring to additional considerations but 
proceeding only from the given systems of PDEs. A systematic computational method 
for constructing the group of symmetries for given system of PDEs has been extensively 
developed by many authors [3]-[7]. A broad review of recent developments in this subject 
can be found in several books as [6], [7]. The methodological approach adopted here is 
based on the symmetry reduction method for PDEs invariant under a Lie group Q of point 
transformations. This means that the groups under considerations are the connected local 
transformation groups acting on the space of p independent and q dependent variables 
E X U, i.e., 



X 



Ag (x, u), u = rig (x, u) ; X = (x^, . . . , x^), u = {u^, . . . , u'^) . (3.1) 



By a symmetry group of a system A of differential equations we understand a local Lie 
group Q transforming both the independent and dependent variables of the considered 
system of equations in such a way that Q transforms given solutions u (x) of A to the new 
ones u(x) of A. The Lie algebra £ of ^ is realized by the vector fields X which for the 
system composed of m PDEs of order k 

A' (x,u('=)) = 0, 1 = 1, ...,m (3.2) 

where u^'^-' denotes all partial derivatives of u up to order k, can be expressed as follows 

X = e^(x,u)9,. +0Mx,u)9„., (3.3) 

where and (f)^ are assumed to be functions of (x, u) & E x U only. The functions ^'^ 
and (j)^ are defined by the infinitesimal invariance criterion [5J 

PT^''UX)A^ =0, l,n = l,...,m (3.4) 

A"=0 

where pr*^'^)(X) is the k prolongation of the vector field. There exist standard algorithms 
for determining the symmetry algebra C and classifying subalgebras Cp of C 
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The symmetry algebra of the MHD equations fl2.1l) -( l23|) . denoted by Ai, has been 
found by Fuchs [8j, and independently by Grundland and Lalague [9]. It is spanned by 
the following 13 Galilean-similitude (GS) infinitesimal generators: 

Pf, = <9x^ , Jk = (^kij [xid^^ + Vjdy^ + Bi&Q^) , Ki = td^^ + d^, , (3.5) 
F = tdt + Xid^^ , G = -tdt - 2pdp + Vidy, , H = 2pdp + 2pdp + Bi&Q^ , 

where ekij is the Levi-Civita symbol, i,j,k = 1,2,3; n = 0,1,2,3. Thus, the MHD 
equations fl2.1l) - fl2.5p are invariant under time (Po) and spatial (Pj) translations, rotations 
(Jfc), Galilei transformations (Ki), dilations (P, G, and H which is the center of Ai; H 
commutes with all the generators). The basis of A4 and the commutator tables are given 
in [9]. In contrast to the results obtained for the (1 + 1) and (2 + 1) dimensional versions of 
the MHD model |8j, the dimension of the Lie algebra A4 for MHD equations fl2.ip - fl2.5l) in 
the full (3 + 1) dimensions (x, y, z and time t) is independent of the value of the adiabatic 
exponent 7 that we treat here as a fixed parameter: 7 > 1. 

A classification by conjugacy classes of r-dimensional subalgebras (1 < t < 4) corre- 
sponding to the GS algebra M. has been established by Grundland and Lalague |9] . Their 
investigation is limited to subgroups of dimension not greater than the number of inde- 
pendent variables. This constraint is imposed by the SRM itself, in order to reduce the 
initial system of PDEs to a lower dimensional system of PDEs which now involves p — r 
independent variables. The method of classification which they have used was developed 
by Patera et al [lOj, and published in a series of papers in the late 1970's. More recently, 
this subject has been treated in such books as [B], [7]. 

The aim of this work [llj was to fulfill the last steps of the SRM applied to the 
MHD equations fl2.1l) - fl2.5p . Obtaining the invariant solutions, called G-invariant solu- 
tions (GIS), under some subgroup Qp C Q requires the knowledge of the group invariant 
I of Qp. For each of the 104 three-dimensional subalgebras of Ai taken from [9], we have 
calculated a set of functionally independent invariants. Suppose {vi,V2,V3} is a basis 
of infinitesimal generators of the Lie algebra £3, then X is a ^s-invariant if and only if 
fj(X) = 0, i = 1,2,3. Thus, for each algebra, we obtain nine functionally invariants. We 
can calculate invariant solutions if the invariants are of the form 

T= {s(x); /i(x,u)}, i = l,...,8 (3.6) 

such that 

In our investigation, we find that the invariants of twelve three-dimensional algebras 
among those classified by Grundland and Lalague [9] do not have the last properties. The 
solutions corresponding to these algebras are called partially invariant solutions (PIS) and 
therefore the SRM is inapplicable. So, we concentrate here only on invariant solutions. 

From the expression given in fl3.6l) . we form the equation 

/,(x,u) = J-i(s(x)), 2 = 1, ...,8 (3.8) 

where J-'i are arbitrary functions, giving the orbits of an invariant solution. From (13.70 
and the implicit functions theorem, we obtain 

u = (s),x) , i,j = l, ...,8 (3.9) 
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where 0j are arbitrary functions. Substituting each of the expressions fl3.9p into equations 
f l2.ip - fl2.5l) . we can reduce the MHD system to a new system composed of ODEs. Such 
system is called the reduced system. The variable s is the symmetry variable and plays the 
role of the independent variable of the reduced system. Solving the reduced system and 
inserting these solutions into (13. 9p . we obtain finally the GIS of the system fl2.ip - fl2.5p . 

In the two next sections, we present ten examples of GIS that we have obtained from 
the GS subalgebras of dimension three [2] (Gj refers to the solution corresponding to this 
algebra given in the title of each subsection). These solutions are also classified by some 
of their physical characteristics and the types of the waves that they may generate. 



4 B=(Bi,B2,0) 

Gi = {Js + Ks + aH, P^, P2}, aeR 

From this algebra, we have obtained the corresponding G-invariant solution 



p = exp 



2a- 

t. 



vi = Uo sin ( — In 
la 



R{t), p = Aj^-^^p, 

V2 = Uo cos 



r ^0 1 




tR{t)_ 





2a 



In 



tR{t) 



Bi = Xo-i / - sin 1 — In ,, 
^' t \2a \_tR{t)\ t)' 



B. = X„,/^cos(^ln 



z 
1 



V3 



(4.1) 

z Wit) 



t 



t 



B3 = 0, 



_tR{t)\ t 

where Uo, Xo and 9o G M. Note that throughout this paper, we denoted the arbitrary 
constants by ( )o- The functions R and W are given by: 

' i?„t(^°'^°-^)exp 



R{t) 



-exp 



Rn 

Ro 

-exp 



2a 



■ 2«2 



^ Wo-—{2Ao + Xl){l + Ht]) 
-Vr„--.Y„^(l + luW) + p— 



for 7 = 1, 
for 7 = 2, 

^(1-7) 



for77^1,2: 



W{t) 



Wo + a {2Ao + Xl) ln[t] for 7 = 2, 

Wo + 7^^t(2-7) + aXl ln[t] for 7 ^ 2 , 
(2-7) 



with Ro e M+, a G M/{0} ; < < l/Aa^ for 7 = 1 and Ao G M+ for 7 > 1, 
Wo G M with sgn[Wo] = sgn[a]. Under these conditions and for t > 0, the solution Gi is 
non-singular and tends to zero for large values of t. 



This solution can be interpreted physically as a nonstationary and compressible fluid 
for which the shape of the flow is a helix of radius Uo and pitch angle = arctan [vs/t/o] . 
Indeed, the motion described by the solution Gi is a circular motion in the x-y plane 
(where lies the magnetic field B) and an accelerated motion parallel to the 2;-axis resulting 
from the action of the Lorentz force: 



1 



2V(|BH=-aX,^|e3, 



(4.2) 
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The Lorentz force Fm acts perpendicularly to B, causing compressions and expansions of 
the distance between lines of force without changing their direction, as does a magnetoa- 
coustic fast wave (F) which propagates perpendicularly to B in an ideal fluid [12] . There 
is no tension force acting on the magnetic fields lines. So, the classical Bernoulli theorem 
is still valid with the pressure of the fluid is replaced by p + (|Bp/2) which is then the 
total pressure. The solution Gi has vorticity that lies in the x-y plane, and referring to 
f l2.1ip . the vortex lines move with the fluid. Consequently, by virtue of Kelvin's theorem 
the circulation 



V d/ 



(4.3) 



c 



around a fluid element is preserved [13J. Also, the solution Gi is characterized by the fact 
that (B ■ V)v = 0, which shows that the velocity of the fluid v is constant along each 
line of force but also meaning that B/p is conserved along the flow, and consequently the 
magnetic field lines remains inextensible [c.f., fl2.10p ]. The current density induced by the 
magnetic field B is equal to 

J = -i [aB2 + Bi] ei + ^ [aBi - B2] es , (4.4) 

which lies in the x-y plane and always has a nonzero component. 



With respect to this algebra, we have calculated the corresponding invariant solution 
Ro 



P 



exp 



(1-7) 



t2(t + ai 

vi = y - ^sm9{z,t), 
z 2a2 



{t + ai) 



{t + ai) 



V2 = - - — cos 9{z,t) 



V3 
Bi 



(t + ai) (t + ai) 



W{t) 



-m\t + ai 



sin6'(z,t), B2 



(t + ai) 



2a2{t + ai) ' 
cos^(z,t), B3 = 0; 



(4.5) 



where Ro-, Ao G M^; Wo-, Uo and Xo G M. The expressions for the functions R and W 
depend on specific values of 7. 

For 7 = 3/2: 



m 

Wit) 
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2a2A<,/(ai)3/2 



exp 



A 

^"2 ^ ^LX2^0 



4a2A arctanh 



1 + 



(t + ai) 



larctanh 



1 + — 

Oil 
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For 7 = 2: 



Rit) 



1 



W(t) = ^In 



exp 
1 

ait 



'2 + In [l + f] 



And otherwise for 7 7^ 3/2, 2: 

(1-7)^(3-27) / 

= \n o„N 2^1(3-27,7,4-27,— 



i?(t) 



(3 - 27) 

* W{s) 
(s + ai) 



a 



-7^(4-27) 



a J (4 - 27) 



Fi (4 -27, 7, 5-27,— 



-ds 



where 2-^1 denotes the hypergeometric function of the second kind. The expression for 
the function 6 is given by 



2ao 



2a2{t + ai) 



(t + ai) 



0oe 



(4.6) 



Solution G2 is always real and non-singular if t > 0, ai > and 02 7^ 0, and tends 
asymptotically to zero when t — » 00. It describes a nonstationary and compressible flow 
in (3 + 1) dimensions with vorticity that lies in the x-y plane. The current density induced 
by (14.51) is equal to 

[a2B2 + Bi]^ , [«2Bi-B2]^ ^^^^ 



ei 



and the resulting Lorentz force takes the form 



(t + ai) 



-62 



(t + ai) 



■peg . 



(4.^ 



Since Fm is a conservative force (i.e., it can be derived from the gradient of the magnetic 
pressure |Bp/2) implies, by virtue of Kelvin's theorem, that the circulation of the fluid 
is conserved. 



Gs = { J3 + P3 + ct^G + a^H, K^, K2], ai 7^ 0, ^2 G 



From this algebra, we have computed the corresponding analytical solution 



W 

X Uo . 



exp 



2(^-1) ^(s) 



P 



A 

^-2a2/ai 



exp 



2(^-7)^(^) 



vi = — sm |t/(s 

Bi = ^t-°2/"iexp 

B2 = ^t-"2/°iexp 
W 



y Uo 

V2 = COS 

t t 



/«2 _ n 

V"i "1/ 



^(s) 



[e{s)l 

sin [^(s)] , 



cos 



1 



y3 = -[W- 1/ai] 



B, = , 



(4.9) 
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where Ro, Ao, Uo, Vo and Xo are arbitrary constants. The functions T and 9 are given by 
^= 9 = 9o + — --, 9,eR, (4.10) 

which both depend on the symmetry variable s = z + (lnt)/ai. The unknown function 
W is determined by solving the reduced form of the equation (12. 2p along the z-axis: 



as Oil Ro 

Rn 



7-— + 2 7-^ 
as \ «! 



Vr(2-7) exp[2(2-7)J^(s 



+ 



^2 



exp 



4--) Hs) 



ds \ai ai 

Equation (14. lip is difficult to solve completely. However, with the ansatz 

W = Cis + C2, Ci, CsGM, Ci^O, 
we can transform (14. lip into a polynomial equation 



ai Ro 



(Ci + 2)7 - 



2^2 



V2 

y^m o 

Ro 



1 a2 



0. (4.11) 



(4.12) 



where m = (2 — 7)(Ci + 2)/Ci and n = 2(2ai — 1)/ (aiCi). By equating the coefficients 
of the same power of W, we obtain five cases of invariant solutions of Eqs (I2.ip -( l23|) . 



Case 1.) Ci = 1, 7 = 4/3, 02 = "i + 1 
p = Rot-'/-^ -3"i)/"i , p = 

2(ai — 1 

y Uo 



,^-2(ai+l)/ai py2{l-ai)/Q!i -j^g^ 



vi = sm t^i, 

t t 



V2 = ---COS^„ V3 = - 



Wi - 



Bi =Xot-("^+^)/°^sin^i, B2 = Xor(°^+^)/"^cos^i, 63 = 0, 
where Ro E Uo, Xo G M; ai > 0. The functions Wi and 6^1 are given by 



9i = 9o In 



C2, 9o E 



Case 2.) Ci = 1, 7 = 2a2/3, ai = 1 

p = Rot^^'-"''^ Wf'\ P = Aot-^''' 



Vl 

Bi 
Bo 



X _Uo 

t t 

V Rn 



sin 9o 



V2 



t t 



COS fc'2. 



a/(2^ 
V(2-a2) 



V3 



■[W^2-l] 



r"Wi"'"^^cos^2, B3 = 



(4.14) 
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where Ro, Ao G M+, Uo, Xo G M; 3/2 < a2 < 2. The functions 1^2 and take the form 



W2 = ln[t] +Z + C2, 02 = 0o- In 



t 



C2, 0o G 



Case 3.) Ci = 1, 7 = 4/3, ai = 1 

p = (2 - «2) + X3 |^^(2"2-5)^ ^ ^ ^^^-2a, ^^(2"2-4) ^ 

^ y Uo Ir . 

vi = - - — sin6'2, V2 = - - — cos6'2, V3 = -[Py2-lJ, 

Bi = Xot-"2 1^2^"'"^) sin ^2, B2 = Xot-""' W!f^'^'^ cos ^2, B3 = , 

where Ao G M^, f/o, G M ; 1 < ^2 < 2. The functions 1^2 and 62 retain the form given 
in solution (14.141) . 



Case 4.) Ci = (2ai - l)/ai , 7 = (2ai + l)/(4ai - 1) 
P 



^2{oi-a2)/oi p;^(2«2-6ai+l)/(2Qi-l) 
(1 - ai)Ro 20:2/01 Ty(2"2-2ai-l)/(2ai-l) 

[2(a2-«i)-l] ' 

a; f/o . . y Uo 

vi = - - — sm 6'3, V2 = - - — cos e's, V3 



(4.16) 



W3- — 



Bi 
B2 



^ . H/,("2-2oi )/(2oi -1) 

y^{2ai - «2) 



3 , 



A/(2ai - 02) 



/oiTT.(a2-2oi)/(2oi-l) 



^-«2/oip^_ 



COS 513, B3 = 0, 



where i?o G M+, f/o G M. The parameters ai and 02 he in these intervals: 1/2 < ai < 1, 
2ai < a2 < (2ai + l)/2. The functions W3 and 6*3 take the form 



(2ai - 1) 



ln[t] 



+ z 



ai (1 — 2ai) 



Case 5.) Ci = (4ai - 2)/3ai , 7 = 6ai/(5ai - 1) 

p = Jl ^2{ai-02)/ai|y^(3a2-8oi+l)/(2Qi-l) 

_^-2o2 /oi p;j^3{o2 -2oi )/(2oi -1) 



P 



B. 



Ro 



2(2ai - 02) 



X Uo . ^ y Uo „ 1 

Vi = sm6'4, V9 = cosc'a, V3 = - 

t t ' ^ t t ^' t 



Bi 



W4 

ai 



' (ai 


-2)i?o 


(4ai - 


■ 3^2 + 1) 


' («i 


-2)Ro 



74 , 



(4ai - 3^2 + 1) 



^~02/0iy^l 



COS 6'4, 63 = 



(4.17) 
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where Ro G Uo E M. Here ai and 02 lie in: 1/2 < ai < 2 and (4ai + 1)/3 < ^2 < 2ai, 
«! > 2 and ai < ^2 < (4ai + l)/3. The functions W4 and 6*4 are given by 



3ai 



ln[t] 



+ C2, g4 = g, + ^+ ^ . Iniiyj, C2, ^^oG 
ai (2 — 4ai) 



Solutions f l4.13p - P:.17p are real and non-singular if VFj > (z = 1, 2, 3, 4) and t > 0, they 
all tend asymptotically to zero for t ^ 00. These solutions describe a nonstationary and 
compressible flow with vorticity that lies in the x-y plane. Note that for the solution 
fl4.14p there is no gradient of fluid pressure. Except for the solution fl4.13p where Fm = 0, 
applying the expression for the Lorentz force (12.80 we find than that only the first term 
(i.e., the magnetic pressure) yields a nonzero contribution, while the second term (the 
magnetic tension force) vanishes. So, Fm is a conservative force that acts only along the 
2;-axis as explained for solutions G\ and G2- Consequently, by virtue of Kelvin's theorem, 
the circulation of the fluid is constant for solutions (I4.13p - (l4.17p . 



5 B=(Bi,B2,B3) 

G4 = {F + aiG + aa/f, Po, P3}, ai ^ 0, 1, ^2 G M 

By choosing 02 = —1, we were able to solve completely the reduced system and we get 
p = (.), p = A- ^1^1^ csc^ {coe) , 



vi =0, V2 = 0, V3 = Wo (x^ + y^r^^ sin"^ {co9) 
CqYq X , . Yq y 



^1 = 7:;ffr^cot(co^) + — (5.1) 



[x^ + r) {x^ + y^) a/x^ + y2 



where Ao G M^; Co, Wo, G M, with Cq 7^ 0. As a result of the incompressibility 

of the fluid, R is an arbitrary function (but positively defined) of the symmetry variable 
s = x/y. The function 6 is given by 

e = arctan (^^) . (5.2) 

Notice that for a given Cq G N/{0}, the trigonometric functions can be expanded in 
terms of rational functions of x and y. This solution is non-singular ii x,y 7^ 0, and 
y 7^ a; tan [/ctt/cq] , k E Ii. We note that this solution does not depend on the variables z 
and t: the axis of symmetry coincides with the z-axis which is also the direction of the 
stationary flow (the paths of fluid elements are streamlines). Solution (15.10 describes a 
static equilibrium [c.f.,the steady-state form of (12. 2p ]: 

Vp = J X B , (5.3) 

where the current density J = (Ji, J2, J3), given by the Ampere's law, takes the form 

Ji = -|^B3, J2 = -|b3, J3 = esc' {coO) . (5.4) 
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Nothing that J x B = — V(|Bp/2) + (B ■ V)B, so to achieve equihbrium the magnetic 
tension must balance not only the magnetic pressure gradient but that of the fluid pressure 
as well. Also, we have (B • V)v = which indicates that the flow velocity v remains 
constant along B. The flow has vorticity given by a; = (cui, 0^2, 0), with 



[x + CoV cot{co6)] [y - CoXcot{co6). , ^, 

cui = ai — V3, uj2 = ai — V3. (5.5) 

(x^ + y"^) (x^ + y^) 

By virtue of Kelvin's theorem the circulation of the fluid is preserved because dv/dt = 
[c.f., equation (12. lip is identically satisfied]. 

The solution (15. ip can be interpreted physically as the propagation of a stationary 
double entropic wave EiEi (resulting from the nonlinear superposition of two entropic 
simple waves Ei) in an incompressible fluid [12]. In the next two solutions, we present 
such type of wave in the context of cylindrically and spirally magnetic geometry. 



G5 = { J3 + aiG + (X2H, Po, P3}, ai, a2 e M 

From this algebra, we get the following invariant solution in cylindrical coordinates (r, if, z) 

p = exp[2(a2 - 0!i)^]Rir), p = Ao exp [2a2ip - 2a2'&{r)] , (5.6) 
Yr = 0, v^ = 0, v^ = Wo exp [aiip - a2^{r)] , 

Br = —exp[a2ip - a2'&ir)] , 

= ^exp [a2(p - a2-&{r)] tan [Y{r)] , = exp [a2(p - a2'(9(r)] , 

where Ao G M"*", Wo, Xo, Zo ^M.; R is an arbitrary function (with i? > 0) of the symmetry 
variable s = r (its level surfaces represent cylinders), and the function ^ is given by 

r dr' 

d= tan(r)— . (5.7) 

The function y is a solution of the nonlinear first order ODE 

AY 

r-, a2/3or^ cos^iY) - = , (5.8) 

dr 

where the parameter /3o is defined by 

(5.9) 

Equation (15. 8p does not have the Painleve property, and is difficult to integrate. If such 
an equation has the Painleve property, i.e., if its general solution has no "movable" 
singularities other than poles (such as essential singularities or branch points), then it 
can be transformed into one of transcendents forms (Ince [H]) and integrated in terms of 
some known functions. The test verifying whether a given ODE satisfies certain necessary 
condition for having the Painleve property, is algorithmic [15] and can be performed using 
a specifically written MATHEMATICA program p^. A complete analysis of this ODE 
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would take us beyond the scope of this paper. However, we can discuss some interesting 
physical properties of this solution. 

The Ampere's law shows that three components of current, J,-, J,^ and J^, are nonzero: 

J, = — B„ J^ = ^tan(r)B„ J, = ^B^ tan(y) - ^B, . (5.10 

r r cos^(y) dr r r 

For ^2 = 0, the current density J is identically zero. So Vx B = permits the expression 
of the magnetic field as the gradient of a magnetic scalar potential, B = — V0m- 

Notice that the equation (15.81) corresponds to the reduced form of the steady-state of 
equation (12. 2p along both the r-and (/J-axes: 

— [p + |B| 2/2] = (B ■ V). B„ b + |Br/2] = (B ■ V)^ B^ , (5.11) 

which both describe a balance between the pressure gradient Vp and the Lorentz force 
Fm- The two terms on the left side of the equalities represent the fluid pressure and the 
magnetic pressure, while the terms on the right side represent the tension forces generated 
by the curvature of the magnetic field lines. Both the magnetic field B and current density 
J lie on constant pressure surfaces since B ■ = and J ■ = 0. 

Solution G5 is non-singular when Y 7^ {2k + l)vr/2, k E 1^, and r 7^ 0. This solution 
is independent of z and t. The solution G5 represents an equilibrium two-dimensional 
configuration for which the stationary and incompressible flow is along the symmetric 
axis (i.e., the ^-axis) of a cylindrically symmetric magnetic field B. Also, (B ■ V)v = 
confirms that v remains constant along B. The fiow has a vorticity field given by 

w = (—v., — tan(F)v„ o) , (5.12) 
\ r r / 

and the circulation of the fiuid is preserved since dv/dt = 0. 



Finally, if we impose the condition that 

B-Vp = 0, (5.13) 

then we can determine the explicit form of the density: 

p = Roexp [2(^2 - ai)ip - 2a2^(r)] , G M+ . (5.14) 

Equation (15.131) simply states that the lines of force lie on the surfaces p = const which 
coincide with the isobaric surfaces. So, any given smooth equilibrium function B, v, p 
and p in M.^ defines a distribution of magnetic surfaces ip = const in M.^ . 
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Ge = { J3 + + G) + a^H, P^, P3}, «i ^ 0, ^2 G M 



In cylindrical coordinates (r, (^9,2;), the corresponding GIS to this algebra has the form 

p = exp[2{a2 - ai)(p] R{s), p = Aoexp [2a2V5 - 2a2^^2(s)] , (5.15) 

Yr = 0, = 0, = Woexp[ai(p - a2'&2{s) 



exp[a2<^ - Ms) - a2Ms)j .... . 
= ^° cos[F(s)]-a,sin[y(.)] '""'^^^'^ 

R V exp[a2y^-^i(g) -a2^2(g)] . r^. r 7 r w 

cos[y (s)J — ai sin[y (s)J '- 

where Ao G M"*", Wo, Xo, Zq G M; and i? is an arbitrary function (but positively defined) 
of the symmetry variable s = re~"^'^. The level surfaces of the symmetry variables are 
logarithmic spirals. The functions i)i and i)2 are given by 

^ ^ r cos(y) ds' ^ ^ r sin(y) d^- 

^ y [cos(r) -aisin(F)] s' ' ^7 [cos(F) — ai sin(F)l s' 
The functions Y solves the integro-differential equation: 

- 7^r\ - "1 sin(r)]2exp [2^1 (s)] - + = , (5.17) 

ds (1 + aj) (1 + ) 

where (3o retains the form (15.91) . and with the conditions that Y 7^ arctan(l/ai) + kn, 
G Z, and s 7^ 0, which ensure that the solution Gq is non-singular. This solution does not 
depend on the variables z and t. Equation ( 15.17^ corresponds to the same reduced form of 
Vp = J X B along both the r- and yj-axes [c.f., equation (15.111) ]. The determination of the 
general solution of equation (I5.17P is a difficult task. We have not able to find an explicit 
solution of (15.171) and we must look for a solution by numerical methods. Nevertheless, 
we can resume some physical properties of this solution. 

The solution Gq represents a two-dimensional spirally symmetric configuration in static 
equilibrium; the steady fiow is parallel the axis of symmetry (i.e., the 2;-axis), and is 
conserved along the magnetic field lines since (B ■ V)v = 0. The vorticity of the fiow is 
given by LJ = {ur, ^^ip,0), with 

uj - '^1 cos(y)v^ ^ _ ai sin(y)v^ 

r [cos(F) -aisin(y)] ' '^'^ ~ r [cos(F) - ai sin(y)] ' ^ ' ^ 

and the circulation of the fiuid is preserved since we have dv/dt = 0. 

The current induced by B has three nonzero components Jr, and J^, given by: 

J _ C(2 cos(F)B2 J _ sin(F)B^ 

^ r [cos(F) — ai sin(F)] ' r [cos(F) — ai sin(y)] 

1/ N rcos(F) + 0:2 sin(y)l , , 

J, = - B^ - a2Br) - \ \J + a A) 

r' r [cos(F) — «! sm(F)J 

2, r , [cos(y)B, + sin(y)B^] dY 

+ [l + a^) exp[-aiif] — — . -77 • (5.19) 

[cos(r ) — «! sm(r )J ds 



13 



If we impose the condition that B ■ Vp = 0, then we get the expression for p 

p = Roexp [2{a2 - a^)^ - 2«2^2(s)] , i?^ G M+ . (5.20) 

So B, V, p and p define now a distribution of magnetic surfaces tl) = const in M.^. 

It is noteworthy that for a2 = 0, the equation (15.171) can be easily solved to give: 
p = exp[- 2ai(p\R{s), p = Ao; v^ = 0, v<^ = , 
= iy„r°i/(°i+i)exp 



ai 



exp 



exp 



(«? + 1) 

.KTI) 



\oi\ sin(F) — cos(F)] 

cos(y) , 

sin(y) , B, = Z„ 



(5.21) 



where Ao G M^, a\ 7^ 0, Wo-, Xq, Zo G M; and R is an arbitrary function (with i? > 0) of 
s = re~"^'^. The expression for the function Y is given by 



Y = Yo + r^^T ln[r] - 

(a? + 1) ^ ^ (a? + 



1) 



(5.22) 



Here we have an example of a two-dimensional configuration with spirally symmetric 
magnetic geometry. Since J = 0, the magnetic field B is potential. 



Gr= {J3 + P3, Pi, P2} 



The reductions of this algebra lead us to this invariant solution of physical interest: 

p = Po, (5.23) 



P = Po, 

vi = ^ sin [V{t) -z], V2 



U{t) 



cos [V{t) -z], V3 = Wo 



/Po ^Po 

Bi =X(t)sin [F(t) -2], B2 = X(t)cos[F(t)-z], B3 = Z„ , 
where po, Vo G ^\ Wo, G M and 7^ 0. The functions [/, V, X and F have the form 



U{t) 
X{t) 



Eo - ^El - po5l 



sm 



2Zo_ 



Eo + /eI^^M! sin (^0o - 



1/2 



.1/2 



V"(t) = ^^0 + 1^0^ + arctan 

/p^^o 



Y{t) = V + arccos 



UX 



where 0o; Vo and l^o ^ I^- The constant E'o is equal to 

Eo=^-^[U{t)f + \[X{t)f 



(5.24) 
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which corresponds to the sum of kinetic and magnetic energies perpendicular to the z-axis. 
The constant 5o is given by 

5o = viBi + V2B2 , (5.25) 

with Si < Ell po ensuring that the solution G^ is always real. Thus, the longitudinal 
component of v along the magnetic field B is a constant quantity. 



The solution (15.231) describes, along with the relations (15.241) and (15.251) . the main 
properties of the double Alfven-entropic wave AEi which results from the nonlinear su- 
perposition of an Alfven wave A with an entropic wave Ei [12] . Such waves of an arbitrary 
amplitude were found for the first time by Alfven p] as nonstationary solutions of MHD 
equations for an incompressible medium (producing no density or pressure fluctuations). 
The wave motion can be explained by the intrinsic nature of the Lorentz force 



(B ■ V)B = -Z0Q2 ei + ZoBi e2 , 



(5.26) 



which corresponds to tension force along the field lines of B. So that the magnetic field 
lines twist relatively to one another, but do not compress. Namely, the double wave AEi 
modifies the flow direction that is not parallel to the magnetic field B. The flow has 
vorticity that lies in the x-y plane, and according to the Kelvin's theorem, the circulation 
of the fluid is not conserved because of the presence of tension forces originating from Fm. 
Note that the instantaneous power dissipated by the magnetic force is given by 



cos 



'Po 



zt 



(5.27) 



Indeed, the equation (15.241) describes basic oscillation between perpendicular fluid kinetic 
energy and perpendicular "line bending" magnetic energy, i.e., a balance between inertia 
and field lines tension. For the particular case Eo = ^/fhj\5o\] v are B coUinear: there is 
no coupling between flow and B. Finally, note that the conservation law of energy [2] 



d_ 

di 



S ,2 , P 



|B| 



+ V- 



S ,2 , 7P 



+ |B| 



v+ (B ■ v)B 







turns out to be trivially satisfied. 



Gs= {F + G, Ki + P2 + (xPs, Pi + /3P2}, a em, (3 > 



The corresponding invariant solution has the explicit form 

z {Wot- Up) ^ . 



V2 



^° ■ [a{l3x - y) + (1 - Pt)z] , V3 = [a{f3x - y) + {I - pt)z] 



R{t) ' J' " R{t) 

_ [Zpt + aXo] ^Yo_ _ Zo 

' aR{t) ' ' R{ty ' R{t) 
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where A^, Wo G M+; Uo, K, X^, Fo, e M, a G M/{0}. The function is given by 

R{t) = pWo - [Wo + (3Uo + aVo]t + i?o , (5.29) 
where the constant Ro must satisfy 

[Wo + m^ + oVl 

in order to obtain p > 0. This solution represents a nonstationary and compressible 
flow in (3 + 1) dimensions. Note that this solution is singular when t = [z + a{(3x — 
y)]/ (3z which coincides with a stagnation point in M'^ (i.e., where p — > oo, v = 0), and 
tends asymptotically to zero for sufficiently large t. The fluid is force-free since both the 
gradient of fluid pressure and Fm are zero. Force-free conditions are widely applicable in 
astrophysical environments because forces other than electromagnetic are comparatively 
much smaller [18]. The vorticity fleld of the flow has the form uj = {uji,uj2,uj3), with 

[PVo t + aWo- Vol [Ro -Uo + a'/3Wo - aVot] [Wot + a/3Vo - Uo] 

^1 = W7T^ > ^2 = 7^77^ , 003 - 



R{t) ' ' aR{t) ' ^ R{t) 

Consequently, by virtue of Kelvin's theorem, the circulation of the fluid is preserved. 



Gs= {F + Ks + aH, P2, P3 + f3P^}, /3 > 0, ae 



The reductions of this algebra lead us to the corresponding analytical solution: 



Vl = U{s), 

Bi = /3yt"exp 

B3 = y t"exp 



, rds'' 



(5.31) 



V2 = K>, V3 = ln[t] + Vr(s), 

'ds' 



- 1+a 



-(! + «: 



/ 



Y 

B o .a 
2 = -jt 



-(l + a^ 



/ 



y 



where Ro, Ao, Vo, Yg and Zo are arbitrary constants. Also, we have introduced the function 

f = U{s) - l3W{s) - s + f3, (5.32) 

that depends on the symmetry variable s = (3\n[t] + {x — I3z)/t. The unknown functions 
U and W are obtained by solving the overdetermined system which consists of f l5.32l) and 
these two equations: 



dU dW 
/?— + -^ + 1 = 0, 
ds ds 

ds Ro 



'-f- — h 7 + 2a 
ds 



'''exp 



Ro 



— + 1 + a 
ds 



exp 



;i-7) 

'ds' 

y 



y 
0. 



(5.33) 



(5.34) 
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We have not succeeded in solving completely this system. However, with the ansatz 



W = Cis + C2, Ci, C2GM, CiT^O, (5.35) 

we can calculate / from equation fl5.34p by setting Ci = —1/2 and 7 = 3. Next, by 
substituting this result into equations fl5.32p and (15.331) . we determine U and W. Finally, 
we obtain the following invariant solution of Eqs (I2.1l) - (l2.5p : 

P = -^(2a + 1)(1 + /3^)[i;^ + (1 + /3^)Z„2]t^"e^^"+^) , 

P = -^1^ + (1 + P')Z'o]t'- e^^^^^) , (5.36) 

, ^ 2 . r,(2+,^)/4 ,1 , PiPz - x) + m - C,) 

' (1 + /32) L ^2(l + /?2)t ^ (I + /32) ' 

Bi = 13 Zo , B2 = Fo t'" , B3 = ^(2«+i) ^ 

e = C72-fln[t]-^^^, (5.37) 

and P > 0, C2, C3, K), and e M. Solution ( 15.36^ is always real and non singular if 
^ > and t > 0. Note that this solution does not depend on the variable y. In order that 
p and p > 0, we must have —3/4 < a < —1/2. This solution describes a nonstationary 
and compressible flow in (2 + 1) dimensions with vorticity along the axis of symmetry: 



where 



[te] 'e2, (5.3^ 



which is everywhere normal to the plane of flow (the x-z plane). The current density 
induced by B takes the form 

J = -{2a + l)Zo ( 1, - (1 + /32), 1 ) . (5.39) 

Thus, the Lorentz force related to solution fl5.36p is given by 

P P 



;i + /?2)t 



1, 0, -/?) (5.40) 



and produced no tension along the lines of force. Also, we have (B- V)v = which implies 
that B/p is constant along the streamlines, and thus B remains inextensible [c.f., (I2.10p . 
So, along both the x- and 2;-axes the Lorentz force Fm causes compressions and expansions 
of the distance between the lines of force without changing their direction, exactly as does 
the propagation of a double magnetoacoustic wave FF which results from the nonlinear 
superposition of two magnetoacoustic fast waves F that propagate perpendicular to each 
other in an ideal fluid [12]. Since Fm is a conservative force (i.e., it comes from the gradient 
of magnetic pressure |Bp/2), the circulation of the fluid (14. 3 p is conserved. 
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Gio = {G + cxH,K^,K^}, ae 



For a = 2, the reduced system corresponding to this algebra admits this exact solution 



Rr. 



A 



t2 

x-U. 



P 



^l^-^exp 



Vl 
V2 

Bi 



+ 



-1 



t ■ Zo{W-al 
y-Vo , al Yot-^ 



-exp 



t Zo{W- al 



-exp 



Xr,t- 



{W - al 



-exp 



{W-al)\ 



2(2-7) 

' dz' 
{W - al) 
&z' 
- al) 

B2 



'dz' 

w 



W 



V3 



-2 



{W - al 



-exp 



dz' 



{W - al) 



(5.41) 



where Ro, Ao, Uo, Vo, Xo, Yo and Zo G M; -Ro and Zo ^ 0, and al = Z^/ Rq. The unknown 
function W of the symmetry variable s = 2r is determined by solving the integro-differential 
equation which corresponds to reduced form of equation fl2.2p along the 2;-axis: 



dz Ro 



2(2 -7) -7 
{XI + Y'i) W 



dW 

"d7 

dW 

dz 



exp 
- 1 



2(2-7) 



exp 



dz^ 

w ^ 

dz' 



(5.42) 



{W - al) 







Ro {W-aD^ 

Equation (I5.42p is difficult to solve completely. However, once again, with the ansatz 

W = Ciz + C2, Ci, C2GM, Ci^O, 
we can transform fl5.42p into an algebraic equation: 

(Ci - 1) + ^(m + 1)W- - l^i±^(Ci -1){W- alY = , 

iXg IXq 

where m = [4 - Ci + 7(Ci + 2)]/Ci and n = {2 - 3Ci)/Ci. Then: 
Case 1.) Ci = 1, 7 = 4/3 



P 
Vl 

Bi 



t2(^ + C2)' 

(x-[/o) ^ ZqXq 
t t Ro 

Yo 

t2' 



P 



Ao 



(5.43) 



Xo_ 

t2 ' 



+ 
B2 



^^^^ _ (y - K) ^ ^ Yo 

t t Rn 



V3 



{Z + C2) 



B, 



t2 



where i?o 7^ 0, G M^; C2, Uo, Vo, Xg, Yo, Zo E M., Zo 0. To avoid singularities and 
to obtain p > 0, we must have t 7^ and z > — C2. Under these conditions, this solution 
tends asymptotically to zero when t —* 00. It represents a compressible and irrotational 
flow of a force-free fluid. 
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Case 2.) Ci = 2/3, 7 = 5/4 



1 + (X,2 + F,^)] (^^ + \ P = # f + ) , (5.44) 



where Ao G M^; C2, f/o, -^o, ^o, € M, 7^ 0. The parameter takes the form 



which can be interpreted as the square of the Alfven velocity. The solution f l5.44p is real 
and non-singular if 2; > 3(a^ — C2)/2 and t 7^ 0, and tends asymptotically to zero when 
t ^ 00. This solution describes the propagation of a compressional Alfven wave in an 
ideal fluid [12], [19]. This can be illustrated with the expression of the Lorentz force: 

F„. = § i^^z + C,- al^ [X. + e,] - ^^^j^ . (5.45) 

The first term on the right-hand side of this equation represents the action of a tension 
force along the magnetic field lines that acts only in the x-y plane (they twist relatively to 
one another, but do not compress), while the second term corresponds to the gradient of 
magnetic pressure, acting perpendicular to the x-y plane, which causes compressions and 
expansions of the distance between the lines of force without changing their direction. So, 
the solution fl5.44p describes a basic oscillation between kinetic fluid energy (fluid inertia) 
and compressional (field line pressure) plus line bending (field line tension) magnetic 
energy. The vorticity of the flow is given by 



UJ 



/2 ^ 



3^ 



{^z + C2-al^ (-n,X„,0). (5.46) 



However, by virtue of Kelvin's theorem, the circulation of the fluid is not conserved 
because of the presence of tension force originating from the Lorentz force: the term on 
the right-hand of the equation (12.111) is nonzero. 



6 Concluding remarks 



We summarize the main results achieved in this paper. Using the subgroup structure of 
the symmetry group of the MHD equations f l2.ip -( l275l) to compute G-invariant solutions 
(GIS) of this system of PDEs, we have considered three-dimensional subgroups where 
codim[^] = 1 in the space of independent variables in order to obtain reduced systems 
strictly composed of ODEs. The Lie algebras of these groups are representatives of the 
conjugacy classes of the Galilean-similitude (GS) algebra that have been calculated previ- 
ously by Grundland and Lalague By means of a procedure of symmetry reduction of 
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Eqs f l2.1l) - fl2.5p . we have found some novel (at the best of our knowledge) exact solutions 
which demonstrate the efficiently of the SRM. Such solutions have been discussed by the 
point of view of their possible meaning. 

In Section 4, we present solutions which describe a nonstationary and compressible flow 
in the presence of a magnetic fleld B = (Bi, B2, 0). In particular, the Lorentz force (except 
for solution (14.131) where = 0) acts perpendicularly to B, causing compressions and 
expansions of the lines of force without changing their direction, as does a magnetoacoustic 
fast wave F which propagates perpendicularly to B in an ideal fluid. There is no tension 
force, so the circulation of the fluid is conserved, which is also meaning that the vortex 
lines are frozen into the fluid. 

In Section 5, we obtained several types of solutions where B = (Bi, 82,33). For the 
algebras G4, G5 and Gq we flnd stationary solutions for which the fluid is incompressible 
and its flow is along the axis of symmetry in the context of slab, cylindrically and spirally 
geometry, respectively. Solution G-j describes the propagation of a nonstationary Alfven 
wave in an incompressible fluid: Fm is purely a tension force and the circulation of 
the fluid is not conserved. Solution Gg represents a nonstationary and compressive flow 
of a force- free fluid in (3 + 1) dimensions. Solution Gg describes the propagation of a 
double magnetoacoustic wave FF in a nonstationary and compressible fluid: causes 
compressions and expansions of the lines of force along both the x- and 2;-axes. There is 
no tension force, so the circulation of the fluid is conserved. Finally for the algebra Gio, we 
obtain solutions (I5.43P and (I5.44p . The flrst represents a compressible and irrotational flow 
of a force-free fluid. The second one corresponds to the propagation of a "compressional 
Alfven wave" in an ideal fluid. Namely, it combines the action of tension force along 
the lines of force and the compressions and expansions of those resulting from magnetic 
pressure force. Obviously, the circulation of the fluid is not conserved. 

It is suitable to recall than an approach of group-theoretical type to the equations 
for the fluid dynamics, especially to the isentropic compressible fluid model and to the 
MHD equations, is not new and has been carried out by many authors (see, for example, 
[8], [20]- [22]). On the other hand, a full classiflcation of Lie subgroups of the symmetry 
group of MHD equations in (3 + 1) dimensions was never been obtained before [9], and 
therefore allows us for a systematic approach to the task of constructing both invariant 
and partially invariant solutions of Eqs (12. 10 - 02. 5p . The objective of future study is to 
take full advantage of some recently developed alternative versions of SRM, namely the 
method of partially invariant solutions (PIS) and the weak transversality method (WTM). 

The notion of PIS was introduced by Ovsiannikov [3] and it relates to the case when 
only a part of the graph P/ of the solution u = /(x) is G-invariant. This means that a 
solution u = /(x) is a partially invariant solution if the number 



(where QTf is the orbit of the graph P/), called the defect structure of the function /(x) 
with respect to the group Q, satisfles the condition 



where r is the dimension of the orbits of Q in E x U. The construction of this type 
of solution leads, through the reduction of the original system A, to two differential 



S = dim (GTf) — p 



(6.1) 



< 6 < min (r, q) , 



(6.2) 
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systems of equations: for 6 dependent and p independent variables and A/Q for 
{q — 6) dependent and {p + 6 — r) independent variables. There is no longer a one-to-one 
correspondence between PIS of A and the solutions of A/Q. For one solution of A/^ we 
have a family of solutions of A. In other words, this approach generates more solutions, 
covering much larger classes of initial and boundary conditions than the classical method. 
Nevertheless, there have been very few examples of its application. Construction of PIS 
requires a much more complex procedure than the "classical" SRM, but Grundland and 
Lalague [23] developed an effective algorithmic tool for this purpose. 

The future study of MHD equations will focus on the construction of PIS invariant 
under four-dimensional subgroups (and some three-dimensional subgroups as mentioned 
above) with the defect structure 6 = 1. Under these assumptions the original equations 
can be reduced to the systems A/Q of ODEs, and one PDE denoted by A^ for one 
unknown function. It could be mentioned also that the procedure cited above produces 
many new reducible solutions which would be, paradoxically, more difficult to obtain by 
the SRM, since they are related to lower dimensional subgroups than Q, leading to the 
reductions of A to PDEs rather than ODEs. 



The second component of a future project will involve the recently WTM [21]. This 
new method is based on a "group invariant solutions without transversality" approach 
developed by Anderson et al. [25] which was designed to overcome the limitation of 
the SRM resulting from the transversality requirement, and extends the applicability of 
this classical method. The notion of transversality refers to the condition imposed on the 
vector fields generating subalgebras of a given system. Consider a r-dimensional subgroup 
Qo C Q, or its corresponding subalgebra Co C £, generated by vector fields 

= C (x, u)9:,M + 0i (x, u)(9„i , a = l, ...,r (6.3) 

with the matrix of characteristics of the vector fields Xa spanning the subalgebra Co 

gi(x,u«) = |0Ux,u)-e(x,u)^| , a = l, ...,r, j = l, (6.4) 

We require that the subgroup Qo acts regularly and transversally on the manifold M. = 
E X U. It means that for each point (x, u) G the relation 

rank{e(x,u)} = rank{C^(x, u), ^^(x, u)} (6.5) 

holds. When this condition, which we call the "strong transversality" condition, is not 
satisfied for a given subalgebra, in principle the classical SRM cannot be applied, i.e., the 
rank of the Jacobian matrix (13. 7p is not maximal. It still possible, however, that there 
exists a certain domain Aio C Ai for which this condition is fulfilled. To distinguish 
such cases, Grundland et al. [2^ have introduced the notion of "weak transversality" and 
have shown that subalgebras with this property can still be used to construct invariant 
and partially invariant solutions using a specific algorithm developed for this purpose. In 
order to do it, we determine the conditions on u = /(x) under which rank requirement 
fl6.5p is satisfied and solve these conditions to obtain the general form of these functions. 
Next we substitute the obtained expressions into the matrix of characteristics (16.41) and 
require that the condition rankQ = is satisfied. This provides further constraints 
on the functions /(x). Finally, we solve the overdetermined system obtained from the 
system under investigation subjected to the constraints coming from requirements of weak 
transversality. If solutions of this overdetermined system exist then invariant solutions 
can be constructed explicitly. 
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